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1  Introduction 

The  radiation  balance  at  high  latitude*  it  dominated  by- 
presence  or  absence  of  mow  and  ice  which  affects  the 
albedo  in  a  great  deal  Any  oceanic  processes  that  af¬ 
fect  the  snow  and  ice  cover  in  polar  regions  will  cause 
the  change  of  deposition  of  sun’s  energy  over  the  earth’s 
surface,  and  hence  the  environmental  change.  Besides, 
the  ice  cover,  ocean,  and  atmosphere  exchange  water 
mass  and  energy  in  a  very  complicated  way.  Three  ma¬ 
jor  feedback  loops  (both  poeitive  and  negative),  in  which 
the  upper  ocean  is  a  key  component,  poosibly  exist  in 
a  theoretical  air-ice-ocean  coupled  system  (Chu,  1990) 
Thus,  investigation  of  the  polar  upper  ocean  thermal 
structure  from  observation  becomes  urgent  Up  until 
now.  the  effort  is  limited  to  study  some  individual  data 
set,  e  g  ,  AIDJEX,  MIZEX,  GEAREX,  etc  ,  separately 
No  investigation  is  pursued  for  the  Arctic  upper  ocean 
thermal  structure  from  the  all  available  data  sets.  This 
study  deals  with  the  most  complete  temperature  and 
salinity  data  sets  in  the  Arctic  Ocean,  which  is  the  U  S. 
Navy’s  Master  Oceanographic  Data  sets  (MOODS).  To 
get  physical  insights  from  the  data,  we  use  a  feature 
model  (llaeger,  1994)  to  transfer  each  MOODS  Arctic 
temperature  (or  salinity)  profile  into  a  set  of  several  char¬ 
acteristics  sea  surface  temperature  (SST),  mixed  layer 
depth,  thermocline  depth,  temperature  difference  across 
the  thermocline  These  characteristics  reveal  strong  tern 
porol  end  spatial  variabilities. 

2  Master  Oceanographic  Data 
Set  (MOODS) 

The  MOODS  data  set  is  a  compilation  of  all  reported  ob¬ 
servations  over  all  the  world’s  oceans.  Due  to  the  shear 
sire  (more  than  5  million  profiles)  and  enormous  influx  of 
data,  quality  control  is  a  difficult  task  MOODS  contains 
different  types  of  data:  (a)  temperature  only  profiles,  (b) 
both  temperature  and  salinity  profiles,  and  temperature, 
salinity,  (c)  and  sound  speed  profiles  We  use  near  30,000 
temperature  and  salinity  profiles  during  1970-1988  for 
this  study  (Fig.l).  To  get  spatial  variability,  we  analyse 
the  profiles  in  several  different  seas  and  the  surround- 


Figure  1  -  MOODS  temperature  nml  salinity  wtntim 
(1070-1988). 


ing  areas:  central  Norwegian  Sea  (0®  —  fi®  E,  fi3  -  (57®  V 
Beaufort  Sea  (160®  -  120' IT.  70®  -  85®  AT).  Gulf  of  Boll 
nia  (15®  -  30® E. 60®  -  67® N),  Denmark  Strait  (30® 
22®  W. 62®  -  68®  N).  Greenland  Se*  (30®  W  -  60®  E.  60® 
65®).  and  Churhi  Sea  (180®  -  160®  W. 65®  -  73®N) 

3  Upper  Layer  Structure 

Arctic  upper  ocean  experiences  a  strong  seasonal  var 
ation  Deep  mixed  layer  appears  in  winter  and  sba 
low  mixed  layer  shows  up  in  summer  During  winte 
tremendous  buoyancy  is  lost  from  the  ocean  surface  di 
to  the  heat  loss  and  salt  injection  caused  by  the  ice  freer 
Buoyancy  loss  along  with  surface  winds  generate  stron 
turbulent  kinetic  energy  (IKE),  which  is  transpose 
downward,  and  deepens  the  ocean  mixed  layer  (OMI 
by  entraining  the  deeper  water  into  OML.  On  the  olh« 
hand,  during  summer,  ocean  surface  gains  buoyancy  ths 
damps  TKE  and  causes  the  mixed  layer  to  shallow 
the  downward  buoyancy  flux  from  the  atmosphere  t 
ocean  is  strong  enough,  mixed  layer  is  very  shallow  wit 
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Figure  2  -  January  sigma-t  during  1970-1077  (mixed 
layer  type). 

a  strong  pyenoeline  below  the  mixed  layer  To  investi¬ 
gate  the  seasonal  and  interannual  variations  of  the  ther¬ 
mal  structure,  we  plot  1970-84  January  (Fig  2)  and  July 
(Fig  3)  <r,  profiles  in  the  central  Norwegian  Sea. 

Comparing  Fig  3  with  Fig.2,  we  found  that  (a)  the 
seasonal  variation  is  very  strong,  (b)  the  winter  mixed 
layer  depth  is  near  300  m,  and  the  summer  mixed  layer 
depth  is  approximate  10  m,  (c)  interannual  variability  is 
quite  weak 

4  Upper  Ocean  Feature  Model 

Most  profiles  exhibit  mixed  layer,  thermoclme  (or  halo- 
eline,  pyrnocline),  and  deep  layer  To  grasp  the  major 
features  of  the  profiles,  a  feature  model  has  been  devel¬ 
oped  in  the  Naval  Oceanographic  Office  (Haeger,  1994) 
for  diagnosing  upper  ocean  thermal  structure  from  tem¬ 
perature  and  salinity  observational  profiles  The  model 
was  originally  called  the  Gradient  Model  Several  pa¬ 
rameters.  which  represent  major  characteristics  of  each 
temperature  profile,  are  SST,  OML  depth,  thermocline 
depths,  and  thermocline  temperature  difference.  Similar 
sets  of  parameters  exist  for  salinity  and  density  profiles, 
eg  ,  sea  surface  salinity  (SSS),  pyenoeline  density  dif¬ 
ference,  etc.  The  feature  model  transforms  each  profile 


Figure  3  July  sigma-t  profiles  during  1970-77  (strat¬ 
ified  type). 

into  a  set  of  these  parameters  The  feature  model  con¬ 
tains  three  components;  (a)  first  guess  feature  model, 
(b)  high-resolution  profiles  interpolated  from  observa¬ 
tions,  (c)  fitting  of  high-resolution  profiles  to  the  feature 
model. 

4.1  Gradient  space 

If  we  consider  profiles  in  the  gradient  space,  i  e.,  CJT  = 
dT[z)/dz,Gs  —  dS{z)/dz ,  each  profile  ran  he  repre¬ 
sented  by  the  surface  value  (SST  or  SSS)  plus  the  gradi¬ 
ents.  eg.. 

(SST.GrfO.ziJ.Crfri.xj).  .C»r(  *■-!.*«)] 

for  the  temperature  profiles  Here,  n  +  I.  «  the  number 
of  data  points,  and  *,(»  =1,2,  n)  are  the  depths  of 
the  data  points  For  example.  100  lemperature/depth 
points  would  produce  99  gradient  values  If  the  surface 
value  is  included,  we  have  the  same  amount  of  data  in  th- 
gradient  space  as  in  the  original  data  set  It  is  difficult 
to  easily  glean  physical  meaning  from  large  number  of 
gradient/flection  point  combinations 

4.2  Thermal  Feature  model 

Observation  shows  that  the  ocean  reveals  layered  struc¬ 
ture  (Fig.4).  Based  on  the  continuity  of  T  and  i\T/dz 
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at  interfaces  of  any  two  layers,  a  feature  model  can  he 
constructed  as: 


TEMPERATURE 


GRADIENT 


f<m)(z)  =  G^m)z  +  SST,  *  €  M, ,  0] 


?<"*>(*)  =  fifth  [(4m)  -  Grn))z  +  -  d,4en)] 

+T H-)(-rf1).  re[-d2,-d,) 


f(,A)(*)  =  G{ft{z  +  d2)  +  f<e">(-rf2),  z  e  [-d3,  -d2) 


[(4,h)  -  &t] V + -  ^r)] 

+f(‘h)(-d3),  ze[-dA,-d2] 


It 


fW(z)  =  G(Td'\z  +  «/<)  +  f«r)(-d<),  z  €  [-rfs,  -<U J 


Figure  4  -  Feature  model 


Tf’U)  =  [(4-1)  -  +  W4dl)  -  °] 

+*?'>(-*).  x  e  [-//,-rfs] 

(1) 

where  (f^m>,f<e"),7<,^r<,r\T^T3<'))  «e  modeled 
temperatures  in  the  mixed  layer,  the  entrainment  zone, 
tht  thermocline,  the  transition  zone,  the  first  and  the 
second  deep  layers.  II  is  the  water  depth,  d\  the  mixed 
layer  depth,  d7  the  depth  of  the  thermocline  top,  rf3 
the  depth  of  the  thermocline  bottom,  dA  the  depth  of 
the  top  of  the  first  deep  layer,  and  d5  is  the  bottom  of 
the  first  deep  layer  (Fig. 4).  Here,  we  assume  constant 
vertical  temperature  gradients  in  the  ocean  mixed  layer 
(Gy"'  s;  0),  in  the  thermocline  (G^),  and  in  the  first 
deep  layer  (G^11*)  and  linearly  varying  with  z  in  the  en¬ 
trainment  zone,  the  transition  zone,  and  the  second  deep 
layer  with  average  values  (G^’.Gy^.G^2*).  By  forcing 
this  feature  model  (1)  to  each  observed  profile,  we  should 
have  a  first  guess  of  the  five  depths  (d, ,  d3,  d3,  dA,  ds) 
and  a  high  resolution  of  temperature/depth  points  in 
the  vertical  in  order  to  obtain  the  six  temperature  gra¬ 
dients  (Gj"\ G(;h\ Such  a 
treatment  provides  the  most  important  features  from  the 
observational  data. 


The  NAVOCBANO’s  DBDB5  data  set  is  used  lo  ob¬ 
tain  water  depth  H.  If  the  five  depths  (rf, ,  rf2,  d3,  dA,  dr.) 
are  known,  we  can  divide  the  data  set  HP  =  (zj,l)) 
into  six  parts  (mixed  layer,  entrainment  zone,  thermo- 
clinc,  transition  zone,  first  deep  layer,  and  second  deep 
layer).  For  each  layer  we  fit  (7),  Zj)  to  the  feature 
model  (1),  and  obtain  a  set  of  temperature  gradients 

(G^\G^\G^\^r\^v),C{r7)) 


4.4  Modeled  profile  (MP)  obtained  by 

THE  ITERATION  METHOD 
A  modeled  profile  (A-f  P)  with  0.5m  resolution  can  be  es¬ 
tablished  by  using  the  feature  model  (l)  if  the  five  depths 
{d],d2,d3ldA,d^)  are  given.  In  reality,  these  depths  are 
not  known  prior  to  processing  the  data  and  vary  from 
one  profile  to  the  other.  We  use  the  iteration  method  to 
obtain  the  optimal  M P, 

First,  we  start  with  a  set  of  first  guess  values  of  the 
depths  and  the  three  constant  gradients, 


D<°  )  =  (/?', $',4°’. 

gV°>  =  (GiC"’,  C^'LcV'1'’) 

In  this  study,  we  choose 


(2) 


4.3  HTGH-RESOLUTION  PROFILES  (HP)  IN¬ 
TERPOLATED  FROM  OBSERVATIONS 
Each  MOODS  profile  is  linearly  interpolated  to  high  res¬ 
olution  (Az  =  0.5m),  denoted  7)  =  T(zj)  for  temper¬ 
ature  as  an  example,  where  Zj  =  Zj_x  —  0.5m(2O  =  0). 


/5<0>  =  (20m,  24m,  Tim.  38m,//) 

G(t]  =  (0,  — 0.5°C/m,  — 0.005°G/m) 

For  each  profile  II P,  we  fit  the  feature  mode  (1)  to  the 
data  (zj,  7j)  since  we  know  the  depths  and  the  constant 
gradients,  and  obtain  a  0-th  order  modeled  profile,  called 
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M />,n !  =  (tjiTj0*)  The  root  mean  square  (RMS)  error 
for  mismatch  of  HP  and  Mi3*0*  is  computed  by 


=  (3) 

where  k  =  0  We  expect  RMS*0^  to  be  large 

Second,  we  use  the  iteration  method  to  obtain  opti¬ 
mal  M  P  for  each  H  P  profile.  Each  depth  can  only  be 
adjusted  one  vertical  grid  (Ar  or  -  Ax)  for  one  iteration 
From  the  kth  order  ( k  starting  from  0,  the  first  guess) 
set  of  depths,  D**’.  we  have  242(=  3s  - 1)  different  com¬ 
binations  of  the  depth  adjustment. 


HA  IM  ••••*«« 


Figure  5  -  Temporal  variation  of  central  Norwegian 
Sea  mixed  layer  depth. 


D<m*+,>=D“>  +  <D£>,  (4) 

where 

iop>»(dM+A,f4 k\4y<y;\ 

Use  ( 1 )  to  obtain  242  modeled  profiles,  among  which  we 
pick  up  a  profile  with  minimum  RMS  error  as  the  (k-f  1) 
th  order  set  of  depths,  £>ft+lV  We  repeat  this  procedure 
until  the  minimum  RMS  error  is  achieved.  When  the 
number  of  iterations  exceeds  a  prescribed  number  (400), 
the  process  is  halted  and  the  RMS  error  is  compared 
to  a  criterion  (0.5*C).  If  RMS  <  0.5 *C,  we  obtain  an 
optimal  set  of  depths  and  gradients  If  RMS  >  0  S*C, 
wr  reject  the  feature  model  (1),  i.e.,  the  HP  profile  can 
not  be  fitted  by  the  feature  model, 

5  Mixed  Layer  Features 

1  he  feature  model  (1)  transforms  any  profile  (if  not  re¬ 
jected)  into  a  set  of  five  depths  (di.rfj.dj.dt.dt)  and  six 

gradients  These 

parameters  represent  most  important  physical  features 
and  are  insensitive  to  the  first  guess  values  The  first 
guess  values  determine  the  number  of  iterations  A  good 
selection  of  first  guess  values  reduce  the  number  of  iter¬ 
ations  Fig  5  shows  both  seasonal  and  interannual  vari¬ 
ations  of  the  mixed  layer  depth  («/t)  in  the  central  Nor¬ 
wegian  Sea  (64  -  67*  N,  0  -  6*  E).  The  ocean  mixed  layer 
is  usually  deep  during  winter  (150  -  300m)  and  shallow 
during  summer  (<  50rn).  The  horizontal  distribution  of 
di  in  the  Arctic  Ocean  (for  the  same  time  period)  shows 
large  spatial  variation. 


6  Conclusion 

The  feature  model  depicted  in  this  paper  is  based  on 
both  the  thermal  structure  and  the  statistics  and  demon¬ 
strates  a  good  capability  to  obtain  physical  insights  from 
observed  temperature  and  salinity  profiles  The  modrl 
uses  iteration  method  to  obtnin  optimal  set.  of  the  param 
eters;  depths  and  gradients  The  model  clearly  shows 
large  thermal  variabilities  and  two-type  structure  This 
feature  model  can  also  he  used  to  prorrss  output  from 
any  dynamical  model  and  to  obtain  meaningful  prod¬ 
ucts,  i.e.,  OML  depth,  thertnocline  depth,  thermor.line 
strength,  and  deep  layer  stratification 
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